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Abstract 

An effective low energy action for Yang- Mills theories is proposed, which 
invokes an additional auxiliary field H^ u for the field strength F^. For a 
particular relation between the parameters of this action a gluon propagator 
with a l/p 4 behaviour for p 2 — ► in the Landau gauge is obtained. The 
abelian subsector of this action admits a duality transformation, where the 
dual action contains a Goldstone boson ip as the dual of H^, and corresponds 
to an abelian Higgs model in the broken phase describing the condensation 
of magnetic charges. The Wilsonian renormalization group equations for the 
parameters of the original action are integrated in some approximation, and 
we find that the relation among the parameters associated with confinement 
appears as an infrared attractive fixed point. 
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1 Introduction 



Presently we are still lacking a proper field theoretical formulation, let alone a 
quantitative description, of confinement in continuum QCD. An illustrative picture 
is provided by the Mandelstam-t'Hooft dual superconductor mechanism of confine- 
ment [|I|]. Here it is assumed that monopoles (with respect to the U(l) subgroups 
of SU (3) color) condense in the QCD vacuum. Then it is supposed that chromo- 
electric flux tubes confine chromo-electric charges in the same way as magnetic flux 
tubes confine magnetic charges inside a superconductor, where electric charges have 
condensed. 

It is very difficult, however, to formulate this idea in the context of a local 
quantum field theory, and even more to prove, that this is a consequence of the 
infrared dynamics of QCD. There are indications for condensed monopoles in the 
QCD vacuum from lattice QCD 0, but this does not help to find a description of 
monopole condensation in the continuum. 

Monopoles are known to arise as classical solutions (solitons) in field theoretic 
models like the Georgi-Glashow model (5*0(3) Yang- Mills theory with Higgs scalars 
in the adjoint representation) ||. t'Hooft has pointed out that similar structures 
could also appear in pure Yang-Mills theories, since the role of the fundamental 
Higgs scalars could eventually be played by some composite fields U |3[]. In any 
case the monopoles can be understood, in the framework of such models, as defects 
in space-time of U(l) gauge fields which arise once the unitary gauge is chosen 

It is notoriously difficult, however, to describe such defects (or solitons in general) 
in terms of quantum fields || and hence to study quantities like effective potentials. 
Moreover, even if one assumes the presence of fields with magnetic charges, a for- 
mulation of a quantum field theory is far from trivial: a Lagrangian for U(l) gauge 
fields, in the presence of both electrically and magnetically charged fields, is nec- 
essarily non-Poincare-invariant (a fixed vector has to be introduced as in axial 
gauges) and either non-local or requires the introduction of a second gauge potential 

Si- 

These latter problems can be circumvented by the introduction of dual fields: in 
the case of a scalar field ip with magnetic charge, e.g., the dual field would be an 
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antisymmetric tensor H^ v (in d = 4). In terms of H^ v it is possible to formulate 
Lagrangians, which are both Poincare invariant and local || [| . 

Employing, on the other hand, dual gauge potentials, it is at least possible to 
formulate low energy models for QCD, which realize the Mandelstam-t'Hooft dual 



superconductor mechanism of confinement [|K| |TT|. Clearly, these models are entirely 
based on assumptions and do not allow to make contact with the bare Lagrangian 
of QCD. 

The purpose of the present paper is twofold: first, we propose a low energy 
effective action for SU(N) Yang-Mills theories, which invokes an antisymmetric 
tensor field H* v in the adjoint representation of the gauge group. We show that, for 
suitably chosen parameters, this effective action implies confinement in the sense of a 
l/p A infrared singularity of the gluon propagator in the Landau gauge. The abelian 
subsector of the action admits a duality transformation, and the corresponding dual 
action is the one of an abelian Higgs model with a Goldstone boson in the broken 
phase. Since accordingly the dual electric charge has condensed in the vacuum, it 
follows that the original magnetic charge has condensed in the vacuum. 

The role of the field H^ v in this action is the one of a composite field for the 
field strength tensor F® u . This can be made explicit by interpreting the low ener- 
gy effective action as a Wilsonian effective action, which satisfies Wilsonian exact 
renormalization group equations (ERGEs). The ERGEs allow to relate the low e- 
nergy effective action in a continuous manner to a bare Lagrangian or high energy 
action, in which the field appears only without space-time derivatives and can 
be eliminated by its algebraic equations of motion of the form i?" = F^ v . 

The second purpose of the paper is thus to set up the system of ERGEs for the 
effective action, and to solve them within some approximation, which allows to keep 
track of all the important parameters. We find indeed that the ERGE flow leads 
us from the bare Yang-Mills action at high scales (including the auxiliary field H®) 
towards a low energy effective action, whose parameters are such that confinement in 
the above sense (1/p 4 behavior of the gluon propagator and monopole condensation) 
occurs. We have thus available both a model for a low energy effective action, which 
realizes the t'Hooft-Mandelstam mechanism of confinement, and a formalism, which 
allows us to compute the corresponding parameters from the bare Lagrangian. 
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Let us, to close the introduction, repeat the essential ideas behind the Wilso- 
nian ERGE approach [fT^| : the starting point is a partition function, in which an 
infrared cutoff k for all propagators is introduced. From the partition function one 
obtains the generating functionals Qk{J) and its Legendre transform Tk(^p), which 
are the usual generating functionals of connected and one-particle-irreducible Green 
functions respectively, where, however, only internal propagators with p 2 > k 2 con- 
tribute. For both functionals Qk{J) and Tk(<p) one can write down exact but simple 
differential equations with respect to k, the ERGEs. For very large values of k, the 
functionals become simply related to the bare Lagrangian of the corresponding the- 
ory, the knowledge of which can thus serve as boundary condition for the integration 
of the ERGEs at some large value k = A. Integrating the ERGEs down to k — 
provides us with the full physical functionals as the effective action To((p). With 
the integration starting at some large but finite value k = A, r (y?) corresponds to 
the effective action in the presence of an UV cutoff A, and is obtained as a function 
of the finite parameters of the corresponding bare Lagrangian. Actually, in order 
to study the renormalizability of a theory, the A dependence of T ((f) can also be 
studied using ERGEs (which now correspond to exact differential equations with 
respect to A JT3|). In the recent years, much progress has been made in employing 



this approach both for scalar theories and gauge theories [Q fT9|j . 

A particularly useful aspect of this method is the fact that it can be extended 
in a straightforward way to functionals of local composite operators flT], |18| : these 
can be introduced in any theory by introducing sources for the composite operators 
in the partition function. After a Legendre transformation one obtains an effective 
action which depends on both fundamental and composite fields. This effective 
action can also be studied using ERGEs. Its dependence on the composite field(s), 
at the starting scale k — A of the ERGE flow, has to be trivial: the action T\ 
depends on the composite field just through a quadratic term without derivatives, 
which would allow - in principle - for its trivial elimination by its algebraic equation 
of motion. The dependence of the full physical effective action r at k = on the 
composite field depends on the dynamics of the theory and is obtained as before by 
integration of the ERGEs. If the full propagator of the composite field, as obtained 
from the quadratic terms in T , possesses a pole in p 2 with the correct sign of the 



4 



residuum, the composite field corresponds to a propagating bound state. Even if 
such a pole is absent, however, the presence of the composite (auxiliary) field can 
be helpful for simplifying the description of the dynamics of the theory. 

Below we will introduce a composite field H^ v for the field strength operator 
F"(A) in non-abelian Yang-Mills theories. Although we do not expect H* v to 
correspond to a physical propagating bound state in such theories, we will find that 
its presence can be very helpful for simplifying the description of confinement. 

In the next section we will discuss several general properties of an effective action 
invoking a composite field H^ v : its formal definition, and the conditions under which 
one obtains confinement. We will describe the duality transformation (of the abelian 
subsector), which allows to interpret confinement as the condensation of monopoles. 

In section 3 we will study the ERGE flow for the effective action in some appro- 
ximation and find that it interpolates indeed between the bare Yang-Mills action at 
the starting scale k = A and a confining action at low scales. Section 4 is devoted 
to conclusions and an outlook. 

2 Effective Action for Yang-Mills Theories with 
an Antisymmetric Tensor Field 

Let us start this section with the definitions of the generating functionals in 
Euclidean Yang-Mills theory in the presence of an infrared cutoff k, and a source 
K° v coupled to the composite operator 

F* v = d,A a u -d v Al + g tf c A^Al . (2.1) 
We will use the short-hand notations 



J-A = j V p j;{-p) Al{p) etc. ; 

V =J_P . (2.2) 
p (2vr) 4 V ; 

Including the standard covariant gauge fixing, ghosts and the infrared cutoff, the 
expression for the generating functional of connected Green functions Gk(J, X> Xi K) 
is 
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exp (-G k (J, x, X, K)) = J T> reg (A, c, c) exp ( - (S YM + S g + AS k ) 

+J-A + x-c + x-c + K-F + K-K) . (2.3) 



Here Sym is the standard Yang-Mills action 



S 



YM 



- f d A r F a F a 
I J " ^ 1 (a/ V ' 



and S g the gauge fixing and ghost part: 



(2.4) 



S g = J d 4 x A; d v Al + d, c a (5 ac d, + g f a bc 4) c c ] . (2.5) 

The term ASk implements the infrared cutoff for the gauge and ghost fields; it is 
given by 



AS k 



T> r , 



-A«(-p) R%(p 2 ) A a u (p) + c a (-p) R k g (p 2 ) c a (p) 



(2.6) 



The functions and R k modify the gauge and ghost propagators such that modes 



with p 2 <C k 2 are suppressed. Convenient choices are 



RUp 2 ) 



p 2 &nv + 

p -V 2 /k 2 



1 \ \ e~P 2 / k ' 2 

~ ~ ) P» Pv) I _ e _ p 2 /k2 



(2.7) 



The functions R^ip 2 ) and R g (p 2 ) vanish for k 2 = 0, and are finite for p 2 <C k 2 such 
that the full gluon and ghost propagators are infrared finite in this regime. 

K denotes the source for the composite operator and we have used the 



freedom |20[ to add a term quadratic in K 1 which simplifies the boundary conditions 



for the Wilsonian generating functionals Qk and Tk for k — > oo [18 



The index "reg" attached to the path integral measure indicates an ultraviolet 
regularization, which is required to make the path integral well defined. Finally 
the generating functionals will be defined entirely through integrated ERGEs with 
suitable boundary conditions, which corresponds implicitely to a regularization in- 
volving higher derivatives. Its explicit form is, however, not required here. 
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The underlying gauge symmetry constrains the physical generating functionals Go 
and Tq severely via Slavnov- Taylor identities. These imply modified Slavnov- Taylor 
identities for the functionals Q k an d 1^ in the presence of a non-vanishing infrared 
cutoff k [15|, In principle their proper formulation requires the introduction of 
sources coupled to the BRST variations of all the fields in the path integral eq. ( |2.3|) . 
Although we will also make use of these identities (e.g. in order to relate the three- 
gluon, four-gluon and ghost-gluon vertices) we have, for simplicity, not shown these 
sources explicitely in eq. (fO|). In particular they are not required for the ERGEs 
themselves. 

The effective action in the presence of the infrared cutoff k is defined, as usual, 
through the Legendre transform 

T k {A,c,c,H) = g k {J, X ,X,K) + J-A + x-c + x-c + K-H . (2.8) 

Sometimes it is more convenient to subtract the cutoff terms AS k from f k and to 
work with 



F k {A,c,c,H) = T k (A,c,c,H)-AS k . (2.9) 



From the path integral ( |2.3| ) and the Legendre transform ( |2.8j ) it is straighforward 
to derive the ERGE for the cutoff effective action 

d ^ = I I wrA Rk (P 2 h ( V/Ta 1 • (2-10) 

2^ (2tt) 4 \8 <pt(-p) dip m (p)J.. 

Here the fields <fi = (A®, c a , c a , H®) denote all possible field appearing as arguments 
of F k or T k , the index i corresponding to the field type and the Lorentz and gauge 
group indices. 

The matrix has non- vanishing matrix elements only in the subsectors (A*, A") 
and (c a ,c a ), cf. the form of AS k given in eq. (|2.6|). The inverse functional 



(5 2 T k /8ip£5(p m )j^, however, has to be constructed on the complete space spanned 
by (Ap, c, c, H^) including the auxiliary field H®. 

It is clear that, in general, T k contains terms with arbitrary powers of the fields, 
and with arbitrary powers of derivatives or momenta p. Moreover, the modified 



Slavnov- Taylor identities [15, 16| require, in general, the presence of all terms, which 
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are invariant under the global part of the gauge group (and have vanishing ghost 
number). Let us, to start with, write down explicitely the BRST invariant terms in 
T k which a) are quadratic in the field strength F" , the auxiliary field H£ u or the 
ghosts, b) contain the lowest non-trivial number of (covariant) derivatives: 



+ -(D^H^f + ^(D fl H, u ) 2 + Z g d,cD ll c 



1 

2a 



+ -L(d II A fl ) 2 + AS k + --- (2.11; 



Here H" is defined by i?" = \s^ vpa and the covariant derivative D^, acting 
on fields ip a in the adjoint representation of the gauge group, by 

D lt < / ? = dr<ffi+gft t A i ll <ff . (2.12) 
The dots in ( p.ll|) denote terms of higher order in the fields, the field strength or 



derivatives, and terms which are not BRST invariant (like a gluon mass term) but 
fixed in terms of the other ones through the modified Slanov- Taylor identities. 

All parameters Z, n, m, h, j3, Z g and a appearing in (|2.11|) depend on the 
scale k, which is equal to the infrared cutoff in AS k introduced in the path integral 
( |2.3|) and present also T k itself. The corresponding ERGEs can be obtained by 
expanding both sides of eq. ( |2.10| ) to second order in the corresponding fields, and 



to the corresponding order of derivatives. These ERGEs will be derived and studied 
in section 3. 

The boundary conditions for the integration of the ERGEs are imposed at some 
large scale k = A, where we require T k =A to correspond to the bare Lagrangian of the 
theory up to additional terms required by the modified Slavnov- Taylor identities. 
Implicitely this leads to a physical effective action T k=Q , after the integration of the 
ERGEs, which contains all quantum contributions with an ultraviolet cutoff A. (In 
some cases it may be desirable to work with a perturbatively improved action at 
k = A [jnj; this will not, however, be employed here). 

The boundary conditions concerning the dependence of r\ = A on the auxiliary 
field H* are fixed by the way the source if" is introduced in the path integral (|2.3|) 
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18| . The corresponding terms exp(K ■ F + K ■ K) can be rewritten with the help of 
an auxiliary quantum field H'\ 



exp(K ■ F + K ■ K) =M J VH'exp(-~(H' - F) 2 + K ■ H'^j . (2.13) 

Now the auxiliary field appears on the same footing as the fundamental fields A, c 
and c, with an action Sh> given by Sh> = \{H' — F) 2 . (This way of introducing the 
auxiliary field H* is close to the field strength formulation of Yang-Mills theories 
f2~I|| .) It is actually convenient to rescale the auxiliary field by a power of A, the 
only scale in the theory, in order to give it the appropriate dimension of a bosonic 
field in d = 4. Then one finds for the boundary condition of the action IV 

T k=A = S YM + S 9 + AS k + ^(AH-F) 2 , (2.14) 



with Sym, S g and AS k as in eqs. (|2.4j )- ( p.6|) . For the parameters Z, n, m, h, j3 and 



Z g in ( 2.11 ) this implies 



Z(A) = 2 
h(A) = . 



n(A) = A 
/5(A) = , 



m(A) = A 
Z g (A) = 1 



(2.15) 



Note that, for arbitrary parameters in ( 2.11 ), the field H£ v could be eliminated 
from ( |2. lip by its equations of motion. This would lead to a dependence of T k on 
the field strength F" of the form 



T k {F) = {Z eff /A)[F; v ) +... , (2.16) 

where the dots denote terms of higher order in the covariant derivatives (induced 
by the terms ~ h, (3) and with 



J eff 



[2.17) 



In terms of Z e ff the boundary conditions ( |2.15|) become 



Zeff(A) = 1 



(2.18) 
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as they should. 

Turning back to the case of general parameters in fl2.ll ), and leaving aside 



the gluon mass for the moment, the terms shown in (EOT) allow to obtain the 



full propagators for all fields. The term proportional to h and quadratic in if" 



corresponds to a kinetic term for H^ v . The expression (dpH^ u ) 2 is actually invariant 
under a gauge symmetry of the form 5H® V = d^A*-, and would not be invertible, but 
the additional term proportional to (3 serves as a "gauge fixing term" and allows for 
a finite H£ v propagator (even for m, n = 0). 

The gauge fixing parameter a in ( |2.11| ) can actually be chosen at will, and 
throughout the paper we will work in the Landau gauge a = 0, which is known to 
be a fixed point of the ERGEs Jl9| . 



Special attention in deriving the propagators has to be paid to the term ■ H, 
which induces a mixing between the gluon and the auxiliary field H*\ this 
modifies the gluon propagator considerably. Omitting the infrared cutoff in ASk for 
simplicity, one obtains 



6(fi(-p) 5ipj(p) J Aa Ab V p 



fj, i J *-v 



5 2 r 



5<pi(-p) Stpj[p) t 



-iSab (p P S^ - Pa 5 W ) ■ Pah(p ) 



Aa lib 



PAH ^ = p 2 (Zm 2 -n 2 ) + Z^ ' ( 2 " 19b > 



5 2 T V 1 



, . I = S ab {8 pK 5 a \- 5 p x 5 aK ) P H H,i{p 2 

\5<pi{-p) 8(pj[p)) Ha Hb 

+8ab(Sp K Pa P\- $p\ Pa Pk + &aX P p Pk ~ &aK P p Px) Ph H ^ 

Phh,i(p 2 ) 



hp 2 + m 2 



2 Z(h-f3)p 2 + n 2 

HH ^ P ' ~ {hp 2 + m 2 ){p 2 {Zm 2 - n 2 ) + Z(3p±) ' (2 ' 19c) 



5 2 T 



6<pi(-p) $<Pj(p), 



Sab P 9 (P 2 ) 
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P 9 (P 2 ) = ^ • (2.19d) 



From eqs. (2.19) and ( 2.15|) one easily obtains simple expressions for the propagators 
at the starting point k = A. Let us now have a closer look at the gluon propagator 
Pa{p 2 ) in eq. (2.19a) (which would remain unchanged, if we would eliminate if" 
by its equation of motion). It depends on four parameters Z, m, n and j3 which, 
in turn, depend on the scale k. Let us now assume , that at some scale k c (possibly 
with k c = 0) the parameters Z, m and n satisfy the relation 

Z i k D - = W*c) = • (2.20) 

Then the gluon propagator Pa becomes 

*<fl-*±£ , A^f • (2,!, 

Thus we find that, for p 2 <C A 2 , the gluon propagator behaves like l/p 4 . Although 
the gluon propagator itself is a gauge dependent quantity, it is widely believed 
that such a behaviour is a signal of confinement: G. West [%| has shown that 



a 1/p 4 behaviour of the gluon propagator in any gauge leads to an area law of 
the Wilson loop, and corresponding results have been obtained in the context of 
Schwinger-Dyson equations for Yang-Mills theories ||23|| . Subsequently we will adopt 
the conventional manner of speaking and call eqs. (|2.20|) and ( p. 21 ) a confining 
behaviour. 

Although we have not shown, at present, that such a confining behaviour actu- 
ally appears for some value of k c (see section 3 below), we will proceed with the 
interpretation of its consequences. We wish to show that an effective action Tk c , 
which exhibits confinement in the sense of eqs. (|2.20|) and fl2.21|) , describes a phys- 



ical situation in which monopole condensation has occurred. To this end we wish 
to perform a duality transformation of the abelian subsector of T^. (Note that this 
does not imply that we neglect the non-abelian contributions to the ERGEs). For 
simplicity we will take only one U(l) subgroup into account, the generalization to 
several U(l) subgroups being straighfroward. 
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Neglecting S g , ASk and the terms not shown explicitely in (|2.3|) , and eliminating 
m 2 by the relation (|2.20|) , the abelian projection of Tk c becomes 



T ha = ±(ZFr-nHr) 2 + ~(dpH) 2 + £(d ll H ia ,) 2 ■ (2.22) 

Here Fp V is the abelian field strength, and all gauge group indices have disappeared. 
For the duality transformation we will also omit the "gauge fixing" term ~ (3, since 
its presence would complicate the duality transformation considerably 

The equations of motion for the abelian gauge field and H^, respectively, 
are then of the form 

- — ' Ti 

h s^vpc d p d\ H Xa — F^ + -gHpy = . (2.23) 
The Bianci identities corresponding to the fields Ap and are given by 



dpd u Hp V = . (2.24) 

Now we introduce dual fields: the dual of A^ is an abelian gauge field Bp with field 
strength F^ u , and the dual of the antisymmetric field Hp U is a scalar ip (with a "field 
strength" dptp). The duality transformation mixes the original fields Ap/Hp V and 
the dual fields Bp/(p in a non-trivial way. It is given by 



d u tp — mBy 



c^-pvpcr \Z Fp a TL Hpu^j , 

Vh dp Hp V 



(2.25) 



with m = n/2\J Zh. 

In terms of the dual fields the action IV becomes 



1 /„ M 2 1 ^ o n2 



T kc (B,<p) = -{F») +-(dp V -mBp 



(2.26) 



It generates the equations of motion for Bp and ip, respectively, 
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^F*+m(d u <p-mB v ) = , 

D(p - m dp Bp = . (2.27) 



After inserting the duality transformations ( |2.25| ) one finds that their content is 
equal to the Bianci identities ( |2.24j ). The Bianci identifies corresponding to the dual 
fields Bp and cp are given by 



d,F* = , 
Spupadpd^ = . (2.28) 



Again, after inserting the duality transformations (|2.25 ), one finds that their content 
is equal to the original equations of motion ( p.23| ). 

Let us now interpret the dual action T^, Q2.26Q : the scalar field <p couples like 
a Goldstone Boson to the dual ("magnetic") gauge field Bp, i.e. we are describing 
the Goldstone and gauge field sector of an abelian Higgs model, where the dual 
"magnetic" £7(1) is spontaneously broken. The dual action Tk c can be obtained 
from a full abelian Higgs model, 

?\ = \ (F*) 2 + (dp + ieBp) 0* (dp - ieBp) <j> + V(|0| 2 ) (2.29) 
in the limit where the vev of the complex scalar field is frozen, 

0=Ji e */« , (2.30) 

and the identification ev = m is made. Since the charge e of the complex field <fi (its 
"electric" coupling to the dual gauge field Bp) is to be interpreted as a magnetic 
charge with respect to the original gauge field Ap, the spontaneous breakdown of 
the dual U(l) gauge symmetry due to dual electric charge condensation is to be 
interpreted as magnetic charge condensation in the original theory (in spite of the 
fact that no non-trivial vacuum appeared in the original theory !). 

Abelian Higgs models as duals of a confining low-energy effective actions of 



Yang-Mills theories have been proposed before []10| , |TTJ. Here, however, we have 
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obtained the Goldstone degree of freedom naturally as the dual of an antisymmetric 
field H^, which was introduced originally for a quite different purpose, namely as 
an auxiliary field for the composite operator F£ v . (During the preparation of this 



paper we became aware of a related proposal in |24|). Also the idea of describing 
magnetically charged scalar fields as duals of an antisymmetric field H^ v is not new 
|8|. U , but this approach was not applied to the dual Meissner effect in @. [| . 

In the following we emphasize that we do not have to assume the validity of 
the confining relation ( |2.20| ), which is required for the duality transformation to be 
possible, in an ad hoc fashion: we have a dynamical scheme at our disposal, the 
Wilsonian ERGEs, which allows to obtain Y^ c - within certain approximations - 
from the bare Yang-Mills Lagrangian. This will be the subject of the next section. 



3 The Wilsonian ERGEs 

In this section we will discuss the Wilsonian ERGEs for the parameters ap- 
pearing (|2.11|) , starting with the functional ERGE Q2.10|) for the cutoff effective 



action. Since all parameters Z, n, m 2 , h, (3 and Z g multiply expressions which 
contain quadratic terms in the fields (~ AA, AH, HH, or cc) the corresponding 
ERGEs can be obtained by expanding both sides of the functional ERGE (|2.10|) 
to second order in the fields, and to the appropriate zeroth, first or second order 
in the derivatives acting on the fields. From the general structure of the ERGEs 



n]-[[TjJ one finds that the right-hand sides of the ERGEs for parameters, which 
correspond to quadratic terms in the fields, involve only coefficients of trilinear and 
quartic terms in the fields. Now we will define the approximation which we will 
employ in the following: we will, on the r.h. sides of the ERGEs, take only those 
trilinear couplings into account, which appear due to the non-abelian structures in 
the terms (F® v ) 2 , F^H^ and d^cD^c in ( |2.11| ). These depend on Z, n, Z g and the 
coupling g (in eq. (|2.12|) and in F£ u ), but the Slavnov- Taylor identities fix, in this 



approximation and in the Landau gauge [[0| g to be 

9 = 9o/Z g (3.1) 
with go (which is equal to the "bare" coupling in T^a) independent of k. If one 
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eliminates if" by its equation of motion, and with an appropriate field redefinition 
of A", one easily finds that a physical running gauge coupling g p h ys (defined at the 
3-gluon, 4-gluon or ghost-gluon vertex) is given by 



g P hys =g/\l Z eff = g / [Zg ■ sj Z eff j (3.2) 

with Z e fj as in eq. ( |2.17| ). 

The approximation consists, in particular, in neglecting the contributions from 
the vertices appearing in the terms (D^H^) 2 and (D^H^) 2 (remember, however, 
that h — [3 — at the starting point) and the contributions from terms not shown 
explicitely in ( |2.11| ). The presence of some such terms in is actually required 



by the modified Slanov- Taylor identities |T^, 16]. (Amongst others, the presence 
of a gluon mass term is necessary. The modified Slavnov- Taylor identities imply, 
however, that this term vanishes in the limit k 2 — »• 0). Note that an approximation 
of the ERGEs of the present type does not imply that violates these identities; 
we rather neglect contributions on the right-hand side of the ERGEs of parameters, 
which themselves are not constrained by the Slavnov- Taylor identities. 

With g determined by ( |3.1| ) we have thus obtained a closed system of ERGEs for 
the 6 parameters Z, n, m 2 , h, (3 and Z g . This is the minimal set which allows a) to 
parametrize the bare action (or the bare Lagrangian), and to obtain the correct 1- 
loop j3 functions, and b) to parametrize a "confining" action where the parameters 



satisfy eq. (|2.20|) . Since actually the quartic vertex without derivatives in (F®) 2 
does not contribute to the ERGEs in this approximation, all 6 ERGEs are of the 
same diagrammatic form as shown in fig. 1. 

The type of external fields {(p a , ipb} and the powers of derivatives <9 M acting on 
them depend on the parameter, whose ERGE one is considering: 

- {<p ai (ft,} = {A, A} and O (d 2 ^j in the case of Z, 

- {(p a , ipb} = {A, if} and O (<9 M ) in the case of n, 

- {<p ai (p b } = {if, if} and O {d°^j in the case of m 2 , 

- {ip a , (p b } = {if, ff } and O (d 2 ^j in the cases of h and j3, 

- {(p a , (fib} = {c, c} and O {d 2 ^j in the case of Z g . 

Note that two different Lorentz index structures are possible in the case {if, if} and 
O (d 2 ), which have to be distinguished in order to separate the contributions to h 
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and (3. The internal fields {ipi,ipj} attached to the insertion of d k R in fig. 1 have 
necessarily to be of the type {A, A} or {c, c}, because only for these fields infrared 
cutoff terms ASk exist. We have to take into account, however, mixed propagators of 
the type {(p e , <£{} = {H, A} or {{p c , ip^} = {H, A} etc., which lead to many different 
contributions even if the external fields {<£>j, <£>,,■} are fixed. 

In order to exhibit the resulting ERGEs we employ the following notations: Pa, 
Pah etc. denote the propagator functions shown in eqs. (2.19), with the infrared 
cutoffs 

-p 2 /fc 2 

*V) = i _ e _ p2/fc2 (3.3) 

due to ASk and the form of the functions Rp V (p 2 ) and R k g {p 2 ) (cf. eq. (|2.7|)) restored. 
This is easily done by replacing Z —* Z + R k and Z g — > Z g + R k in eqs. (2.19). 
Actually the exhibited form of the equations is general and allows for different choices 



of the cutoff function (3.3). (Our choice guarantees that all integrals appearing on 
the r.h. sides of the ERGEs are ultraviolet finite; however, we were not able to 
perform these momentum integrations analytically.) 

Because of the required development in powers of derivatives (or external mo- 
menta) the quantities p 2 ^t etc. appear frequently, and we define for convenience 

BP d 2 P 

p '^ h • p "^wf {3A) 

in all cases Pa, Pah etc. Finally we define = k 2 d/dk 2 , and performed the 
combinatories of the gauge group indices in the case of a SU(N) gauge group. The 
6 ERGEs then become 



B k Z 



Ng_ 

16tt 2 



dp 2 p 4 d k R" Z l ( ^-P% + WIP'a + P\P. 



)/ -7 2 I 31 7,3 



2 TV 



?2 r," 



6 



-nZ{^-P 2 A P AH + \p A Pah P' a 



10 

y 



-L—P 2 p' _L p2 pH 
+— r A * AH + r A * AH 



W{PaPah[~P' ah 



— r a 



AH 



7_ 

12 



1 



Pa Pah Pa 



16 



13 



11 



12 



, 1 7 2 p3 



dkn = x& L°° dp2 pA dk Rk n P " A {l z Pa ~ i n Pah 



d k m 2 



7,2 



Ng 
16tt 2 



/ dp 2 p 4 d k R k n 2 P\ 
Jo 



Nn 2 r°° 3 
d k Z a = ^j^ dp 2 p A d k R k Z 2 P A P g -{P A + P t 



9) i 



d k h 

d k f3 



Ng 2 
16tt 2 

-d k h 



f dp 2 p 2 d k R k n 2 P 2 (lp A + A -P' A + \PI 



(3.5a) 



(3.5b) 



(3.5c) 



(3.5d) 



(3.5e) 



Note that all integrals are trivially UV finite since, with the present choice of R k , 
d k R k decreases exponentially for large p 2 , and IR finiteness is ensured by the presence 
of the IR cutoff terms in the propagators. 

Let us first insert the "starting point values" (|2.1|) for the parameters on the r.h. 
sides of the ERGEs (3.5). Then all integrals can be performed analytically, and we 
should obtain the correct 1-loop (3 functions. We find 



This leads, indeed, to 



8 7- Ng » 31 

o k n = ■ -A , 

16tt 2 4 

o 2 ^0 1 A 2 

dkm = 16V-2 A ' 

16tT 2 12 ' 

8 7- N9 ° 3 f3 6) 



d k Z eff = d k [z =— , 3.7 

n \ m? 16tt 2 6 V ; 
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thus we obtain the correct 1-loop (3 function for g p h ys as defined in eq. (|3.2| ). From 
eqs. ( |3.6|) and ( |3.7|) one concludes that Z, n, m 2 , Z g and Z e ff start to decrease with 
decreasing k 2 , whereas h (and hence /3, cf. ( |3.5ej )) start to increase. 

The complete system of ERGEs (3.5) can be integrated numerically, with the 
choice of g = g p h ys (A 2 ) as the only freedom at the starting point. In fig. 2 
we show a typical result obtained with g = 1.3. We plot directly Z e ff(k 2 ) = 
Z(k 2 )—n 2 (k 2 )/m 2 (k 2 ) and, for convenience, a p h ys {k 2 ) = g 2 /(4irZ 2 (k 2 )Z e ff(k 2 )) ver- 
sus -£n(k 2 /A 2 ). We find that at a small, but finite value of k 2 = k 2 c (= 1, 94- 10~ 5 A 2 
in the present case) both Z g {k 2 ) and Z e ff{k 2 ) vanish, thus a p h ys (k 2 ) runs into a 
Landau singularity and the r.h. sides of the ERGEs (3.5) explode. The parameters 
Z(k 2 ), n(k 2 ) and m 2 (k 2 ), however, remain finite and non- vanishing for k 2 — > k 2 (in 
the present case we have Z(kj) = 0.24, n(k 2 ) = .43 ■ A, m 2 {k 2 c ) = 0, 76 • A 2 ). 

From the previous section, in particular from eq. (2.20) and the subsequent 
discussion, we know that a vanishing of Z e ff(k 2 ) corresponds to a confining form 
of the gluon propagator. In the corresponding expression eq. (|2.21|) a dimensionful 
quantity A 2 = m 2 / [3 appears which could be related to the slope of a linearly rising 
gg-potential obtained by a "dressed" one-gluon exchange (with a gluon propagator 
as in eq. (|2.21|) ) in the non-relativistic limit. As such it should be independent from 
the starting scale A. However, A is the only dimensionful scale in our approach, 
hence all dimensionful quantities are necessarily proportional to A. The physical 
meaning of A is implicit in the choice of the bare coupling go : go = g p h y s 

(A 2 ). Given 

the A-loop (3 function for g p h ys (the 2-loop f3 function is not obtained correctly within 
our approximation), the combined dependence of A 2 on A 2 and g 2 should thus be of 
the form 

A 2 = A 2 e -(i6^/ii 9o 2 ) (i + o (&)) . (3.8) 
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We have calculated the quantity 



A 2 M) rn 2 (k 2 ) ^" yj 

for various choices of go, where the pattern of ERGE flow for all parameters is 
always of the form presented in the case go = 1.3 (but with different results for 
k 2 JA 2 , of course): g = 1.3: C = 1.68 • 10~ 2 ; g = 1.4: C = 1.38 • 10~ 2 ; g = 1.5: 
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C = 1.20 • 1(T 2 ; g = 1.6: C = 1, 09 • 10" 2 ; g = 1.7: C = 1.01 • 10~ 2 . Thus, although 
the corresponding values of A 2 (fc 2 )/A 2 vary over 2 orders of magnitude, C remains 
approximately constant with an accuracy in agreement with eq. ( |3.8| ). (For even 
larger values of go the ratio A 2 /A 2 becomes larger than 1, i.e. the UV cutoff A is 
below the "confinement scale" A c , whereas for smaller values of go the number of 
discrete steps for the integration of the ERGEs (~ A 2 /k 2 ) becomes extremely large 
without any new insight being gained). 

Actually the fact that the "confining" relation Z e ff(k 2 ) = and hence a Landau 
singularity in a p h ys {k 2 ) appears already at a finite value of k 2 = k 2 c is to be considered 
as an artefact of the present approximation to the ERGEs, which does evidently not 
yet capture the entire infrared dynamics of Yang-Mills theories ( possible 
gluon condensate). Focusing on the three parameters Z, n and m 2 , which appear 
in the "confining" relation (|2.20| ), it is very interesting, however, that an analytic 
result can be obtained, which is independent of the other parameters: let us consider 
the r.h. sides of the corresponding ERGEs in the case where the relation Z{k 2 ) — 
n 2 (k 2 )/m 2 (k 2 ) = is fulfilled. In addition we assume that the scale k 2 is already 
far below the "confinement scale" A 2 = m 2 /f3, k 2 <C A 2 . Since only momenta with 
p 2 <k 2 contribute to the integrals of the ERGEs (3.5) (remember the exponential 
damping), this assumption allows to simplify the propagators considerably and to 
obtain analytic results. We find 



k 16vr 2 ' k 2 ' 4 
Ng 2 A 2 n 



d k n 



16vr 2 k 2 4 ' 
Ng 2 A 2 m 2 



d * m -T^-¥-T ■ (3 - 10) 



One easily derives 



i.e. whereas the individual parameters Z, n, and m 2 continue to run in a way 
which depends on the running of g 2 , Z e ff = is a quasi-fixed point of the ERGEs. 
Moreover, in agreement with our numerical analysis, this quasi-fixed point is infrared 
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attractive. Were it not for the singular behaviour of g 2 for small k 2 due to the 
vanishing of Z g , the integrated eq. (|3.11| ) would lead to a smooth vanishing of Z e ff 
for k 2 — > 0. Again, however, we cannot check in how far the quasi-fixed point 
Z e ff = depends on the present approximations to the ERGEs. 



4 Conclusions 

The intentions of the present paper are twofold: first, we proposed an effective 
low energy action for Yang-Mills theories, which involves an additional antisymme- 
tric tensor field H® u introduced as an auxiliary field for the composite local operator 
F° We have seen that, for a particular relation between 3 of the 6 parameters 
in this action, a confining gluon propagator in the sense of a 1/p 4 behaviour for 
p 2 — > is obtained. On the other hand, we were able to relate this behaviour 
to a condensate of monopoles: the presence of the field H^ u in the original action 
generates automatically the presence of a Goldstone boson (p in the dual action (of 
the abelian subsystem), which is thus of the form of an abelian Higgs model in the 
broken phase (with frozen radial excitations). This "dual" role of the auxiliary field 
H® is certainly a particular feature of gauge theories in d = 4 dimensions. An 
amazing and conceptually important feature is the fact that, from the point of view 
of the original action, the confining behaviour of the gluon propagator is not directly 
related to a non-trivial vacuum (it can be speculated, however, that a vanishing of 
Z e ff is a necessary condition for the formation of a gluon condensate). 

Second, we have discussed a dynamical scheme, which allows to compute the 
corresponding low energy effective action from the bare Yang-Mills Lagrangian: the 
Wilsonian ERGEs including auxiliary fields. Since this constitutes, a priori, an 
exact formalism, it allows for various approximations in order to turn it into a 
closed system of differential equations. These approximations can be chosen such 
that they correspond to the physical problem under consideration. Here we did not 
solve the ERGEs for the parameters corresponding to higher dimensional or non- 
BRST-invariant terms in the effective action, whose knowledge is not required for 
the check of the "confining relation" eq. fl2.20|). The approximation consisted then in 



the neglect of the contributions of these parameters to the ERGEs of the 6 important 
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parameters Z, n, m 2 , h, (3 and Z g . We have seen that, within this approximation, 
the relation (|2.20[) constitutes indeed an infrared attractive quasi-fixed point of the 
ERGEs. On the other hand, the fact that the relation ( 2.2(J| ) is assumed already at a 
finite scale k 2 = k 2 , and the physical gauge coupling runs into a Landau singularity, 
is obviously an artifact of the approximation. We have certainly not yet solved the 
problem of finding approximations to the system of ERGEs, which allow to obtain 
numerically reliable results in the infrared. 

However, the present approach may point into a direction where this problem 
may be solvable. We have seen that the duality transformation of the abelian sub- 
system provides us with the action of an abelian Higgs model in the broken phase. 
This model has trivial infrared fixed points in the Wilsonian sense, since it contains 
no massless physical fields which could contribute to the ERGE flow in the limit 
k 2 ~ p 2 — > 0. Duality transformations can also be considered for full non-abelian 



actions, see [25| for first steps in this direction. Here the genuine problem is that one 
obtains actions, which are non-local, non-polynomial in the fields, and/or involve 
additional auxiliary degrees of freedom. On the other hand, all these features are al- 
ready present in Wilsonian effective actions in general, and the Wilsonian functional 
ERGE is capable to cope with then. Let us now speculate that the result of a duality 
transformation of a full non-abelian action, again in the presence of auxiliary fields 
H^ u , leads us again to dual actions which contain no massless physical fields. Then 
we would obtain automatically infrared fixed points for the Wilsonian ERGEs for 
the dual action, which allow to obtain numerically reliable results in the infrared. It 
then remains to reinterpret the infrared effective action of the dual theory in terms 
of the original theory in order to extract its physical consequences. Although we are 
aware of the speculative nature of this approach it is presently under investigation. 



Acknowledgements 

We have profited from numerous discussions with B. Bergerhoff, D. Jungnickel, 
A. Weber and C. Wetterich. 



21 



References 



[1] S. Mandelstam, Phys. Rep. 23C (1976) 245; G. t'Hooft, Nucl. Phys. B190 
(1981) 455. 

[2] Lattice 96, Proceedings; Nucl. Phys. B (Proc. Suppl.) 53 (1997). 

[3] A. Polyakov, JETP Lett. 20 (1974) 194; 
G. t'Hooft, Nucl. Phys. B79 (1974) 276. 

[4] M. Blagojevic, P. Senjanovic, Phys. Rep. 157 (1988) 233. 

[5] A. Di Giacomo, M. Mathur, Phys. Lett. B400 (1997) 129. 

[6] K. Bardakci, S. Samuel, Phys. Rev. D18 (1978) 2849; 

E. Akhmedov, M. Chernodub, M. Polikarpov, M. Zubkov, Phys. Rev. D53 
(1996) 2087. 

[7] D. Zwanziger, Phys. Rev. D3 (1971) 880. 

[8] M. Mathur, H. Sharatchandra, Phys. Rev. Lett. 66 (1991) 3097. 

[9] K. Lee, Phys. Rev. D48 (1993) 2493. 

[10] S. Maedan, Y. Matsubara, T. Suzuki, Prog. Theor. Phys. 84 (1990) 130 and 
refs. therein. 

[11] M. Baker, J. Ball, F. Zachariasen, Phys. Rep. 209 (1991) 73, Int. J. Mod. Phys. 
All (1996) 343; 

M. Baker, J. Ball, N. Brambilla, G. Prosperi, F. Zachariasen, Phys. Rev. D54 
(1996) 2829. 

[12] K. Wilson and I. Kogut, Phys. Rep. 12 (1974) 75; 

F. Wegner, in: Phase Transitions and Critical Phenomena, Vol. 6, eds. C. Domb 
and M. Green (Academic Press, NY 1975). 

[13] S. Weinberg, in: Proceedings of the 1976 International School of Subnuclear 
Physics, Erice, ed. A. Zichichi; 
J. Polchinski, Nucl. Phys. B231 (1984) 269. 

22 



[14] G. Keller, C. Kopper, Phys. Lett. B273 (1991) 323; 

G. Keller, C. Kopper, M. Salmhofer, Helv. Phys. Acta 65 (1997) 32; 

C. Wetterich, Phys. Lett. B301 (1993) 90; 

C. Wetterich, M. Reuter, Phys. Lett. B334 (1994) 412; 

T. Morris, Int. J. Mod. Phys. A9 (1994) 2411; 

M. D'Attanasio, T. Morris, Phys. Lett. B378 (1996) 213; 

For more references see J. Cornelias, Y. Kubyshin, E. Moreno, Nucl. Phys. 

B490 (1997) 653. 

[15] C. Becchi, in: Lectures given at the Summer School in Theoretical Physics, 
Parma 1991-1993, eds. M. Bonini, G. Marchesini and E. Onofri, Parma Univer- 
sity 1993; 

M. Bonini, M. D'Attanasio, G. Marchesini, Nucl. Phys. B418 (1994) 81, Nucl. 
Phys. B421 (1994) 429, Nucl. Phys. B437 (1995) 163, Phys. Lett. B346 (1995) 
87. 

[16] U. Ellwanger, Phys. Lett. B335 (1994) 364. 

[17] U. Ellwanger, Z. Phys. C58 (1993) 619. 

[18] U. Ellwanger, C. Wetterich, Nucl. Phys. B423 (1994) 137. 

[19] U. Ellwanger, M. Hirsch, A. Weber, Z. Phys. C69 (1996) 687 and |Ke]> 



ph/9606468| , to appear in Z. Phys. C; 



U. Ellwanger, |hep-ph/ 9 702309] , to appear in Z. Phys. C 



B. Bergerhoff, C. Wetterich, |hep-ph/ 9708421) . 



[20] R. Haymaker, Riv. Nuovo Cim. 14 (1991), No. 8. 
[21] M. Holperu, Phys. Rev. D16 (1977) 1798; 



M. Quandt, H. Reinhardt, [hep-th/9707185| and refs. therein. 



[22] G. West, Phys. Lett. B115 (1982) 468. 

[23] C. Roberts and A. Williams, Prog, in Part, and Nucl. Physics 33 (1994) 477, 
and references therein. 



[24] K. Kondo, preprint CHIBA-EP-99, |hep-th/9709T09 . 



23 



[25] K. Itabashi, Prog. Theor. Phys. 65 (1981) 1423; 

L. Mizrachi, Annals Phys. 153 (1984) 417 and Phys. Lett. B139 (1984) 374; 
A. Caticha, Phys. Rev. D37 (1988) 2323. 




Figure 1: Diagrammatic form of the ERGEs for the parameters corresponding to 
terms quadratic in the fields in T k . The external fields f a , ipi, have to be chosen to 
be A, H, c, or c as described in the text. The crossed circle denotes an insertion of 
the function d k R k (p 2 ). 
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Figure 2: Result of the numerical integration of the ERGEs (3.5) with go = 1.3. 
Z e ff(k 2 ) and a p hy S (k 2 ) are defined in eq. (|2.17|) and below eq. ( |3.7p , respectively. 
The ERGE flow starts at —£n(k 2 /A 2 ) = and stops due to the Landau singularity 



at -in(k 2 /A 2 



10.85. 
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